Reply to Comment by Borisenko et al. on article 'A de Haas-van Alphen study of the 
Fermi surfaces of superconducting LiFeP and LiFeAs' 
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Recently, Borisenko et al. [l| have posted a Comment 
where they suggest an ahernative interpretation of our 
de Haas-van Alphen (dHvA) measurements on the su- 
perconductor LiFeAs 0. In our original paper j3l we 
concluded that our measurements of the bulk Fermi sur- 
face were not consistent with the surface bands observed 
thus far by ARPES. Borisenko et al. [l| dispute this and 
suggest the two measurements are consistent if some of 
the orbits we observe are due to magnetic breakdown. 
We argue here that this scenario is inconsistent with the 
experimental data and therefore that our original conclu- 
sion stands. 

The Comment by Borisenko et al., [H consists of sev- 
eral claims. First, it is claimed that the technique of 
interpreting the dHvA spectra in terms of the calculated 
band structure, and accounting for the deviation from the 
calculated density functional (DFT) bands by assigning 
small mutual shifts to the bands, fails in case of LiFeAs. 
To prove that, they plot together the unshifted bands 
(our Fig. 2c-left) overlapping them with the measured 
data of Fig. 2c-right. In fact, our Fig. 2c-right shows 
that the shifted bands agree with the experiment within 
the experimental noise. We note that the required shifts 
are surprisingly small, smaller than in any other multi- 
orbital metal, correlated or not, e.g., MgB2 or Sr2Ru04 
d, We also note that, contrary to the claim in the 
Comment [J, although in general the hole and the elec- 
tron bands are allowed to shift independently in the ex- 
perimental fitting procedure, the final result is checked 
for physical consistency by ensuring the total number of 
electrons is conserved. In fact, the very small size of the 
shifts needed in the case of LiFeAs renders this issue ir- 
relevant. 

Their second claim is that in addition to the afore- 
mentioned shifts a huge mass renormalization is needed 
to reproduce the observed cross-section areas, depicted in 
Fig. 2c. This is a misconception. Effective masses have 
no effect on the areas, and no other modifications of the 
calculated bands, besides small shifts, had been made in 
making Fig. 2c-right. The authors of the Comment [ij 
were perhaps confused by the fact that dHvA allows for a 



separate extraction of the effective masses (from the tem- 
perature dependence of the peak amplitudes) , and these 
appear to be very orbital dependent, with one particular 
band showing mass renormalization of 4.9. 

The authors of the Comment appear surprised that 
LiFeP (though not LiFeAs) requires up to 73meV shifts 
to describe the experimental results. However, this is 
in fact a very typical number for any multiband mate- 
rial. We recall that MgB2, which is exceedingly well stud- 
ied by dHvA, ARPES, and many other techniques, and 
lacks magnetism, correlations, and other complications 
of ferropnictides, requires shifts of up to 100 meV. On 
the other hand, the parallel between LiFeAs and LiFeP, 
that these authors try to exploit here, is indeed extremely 
useful and important. The bulk sample quality of LiFeAs 
as of now is not quite as good as that of LiFeP, but the 
fact that in LiFeP all orbits below 3 kT are seen and are 
reproduced exceedingly well by our procedure, leaves no 
room to doubt the applicability of the same procedure for 
LiFeAs. The same procedure has also been used to fit in 
very high detail other ferropnictides such as LaFePO 0], 
SrFe2P2 [si, BaFe2P2 [9] and BaFe2As2 We note in 
passing that the authors appear to have another miscon- 
ception, that a band with a large mass renormalization 
would require a larger shift to achieve the same adjust- 
ment in the orbit areas than an unrenormalized band. In 
fact the opposite is true; while we had to shift the DFT 
band 1 by 60 meV, the renormalized band (m*/m = 2.3 
for this band) would have to be shifted by 26 meV. 

Yet another point of the Comment is that the observed 
near degeneracy of the 2a and 5a orbits is accidental. 
This is indeed true, and we say that explicitly in our 
paper Moreover, DFT being only an approxima- 
tion, the real bands may or may not hold this degen- 
eracy. However, the reader should be aware that the 
LMTO-ASA method to which Borisenko et al. [H re- 
fer is only a less accurate method for solving the same 
DFT equations than the LAPW method used by us. In 
other words, a high accuracy solution of DFT equations 
yields this accidental degeneracy, and a less exact solu- 
tion of the same equations violates it. Of course, a more 
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FIG. 1: Possible orbits in the electron bands of LiFeAs at 
kz = 0.5c*. NBD shows the two non-breakdown orbits. The 
orbits resemble a flower with 4 petals. The breakdown orbit 
BDl has the central part plus 1 petal, BD2 has 2 additional 
petals and BD3 has 3 petals. 



accurate approximation than DFT may also violate it, 
but one cannot address this issue in the way the com- 
ment [H does. We also note here that the authors of Ref. 
[l[ do not state if the structural parameters used in the 
unpublished LMTO calculation they cite are the same 
as in ours. This could be another source of discrepancy 
because the band-structure depends strongly on the As 
internal position. 

The main substance of Ref. [ij is the suggestion that 
the dHvA data might be consistent with the ARPES re- 
sults if magnetic breakdown is considered. In the Com- 
ment, the authors compare our data to their unpublished 
ARPES measurements where they have observed both 
electron Fermi surfaces as a function of kz (in Ref. Q 
only one electron surface and no kz dispersion was re- 
ported). In fact, in the comparison shown in Fig. 1 
of their Comment the authors use only the experimen- 
tal Fermi surface cross-sections at the center and top of 
the zone and then estimate the angle dependence of the 
dHvA frequencies using the result of the LMTO calcula- 
tions at low angle only [6], extrapolating the low angle 
result to higher angle instead of using the actual calcula- 
tions. Such a procedure can hardly be considered to be 
rigorous but nevertheless they suggest that the e dHvA 
orbit (see Fig. [3]) is consistent with a magnetic breakdown 
orbit coming from the electron orbits at top of the zone 
and that the 7 and /3 orbits come from the conventional 
(non-magnetic breakdown) orbits from the inner electron 
sheet (at the top and center of the zone respectively). It 
is suggested that this occurs because the gap between 
the electron sheets (which is induced by the spin-orbit 
interaction) is very small. We had carefully considered 
this scenario and concluded that it is inconsistent with 
the data. 

Magnetic breakdown occurs when the electron has suf- 
ficient energy to tunnel between bands. The probabil- 
ity of it occurring depends exponentially on the field, 
p = cxp{—B()/B), where the breakdown field Bq depends 
on the size of the gap. For the electron orbits in LiFeAs 
(see Fig.[T]), for B <C Bo two orbits should be seen corre- 
sponding to the inner and outer electron bands, whereas 
for B ^ Bo a single orbit corresponding to the ellipse 
in the unfolded zone is expected (BD2(i) in Fig. [T|). In 
the intermediate field limit B ~ Bq both types of orbit 
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FIG. 2: Relative probabilities of performing the orbits shown 
in Fig.[Tl Non breakdown = (1-p)*, BD2 = 2/-hV(l-p)^ 
BDl or BD3 = 4p2(l -p)^ where p = exp(-Bo/-B). 



would be seen along with two further orbits correspond- 
ing to the inner orbit plus one or two or three additional 
sections of the outer orbit (see Fig. [T]). Because the el- 
lipse orbit requires 4 breakdown gaps to be crossed, its 
probability is proportional to p**, whereas the low field 
orbits are proportional to (1 — p)* and the intermediate 
orbits, which require two tunneling events and two reflec- 
tions, are proportional to p^(l — p)^. These probabilities, 
multiplied by their frequency of occurrence, are plotted 
as a function of B/Bq in Fig. [21 

The interpretation suggested in Ref. [l[ requires the ob- 
servation of both breakdown and non-breakdown orbits 
so as shown in Fig. [2] we require i? ~ Bq, and hence all 
three types of breakdown orbit should be observed with 
approximately equal amplitude. In additional, break- 
down should be equally likely to occur at the center and 
the top of the zone and hence the number of breakdown 
orbits is doubled. Using the estimated ARPES F{9) 
curves from Ref. jl| we plot all these predicted orbits 
along with the data in Fig. [3] It can be seen that only 
the breakdown curve BD2Max is consistent with the data. 
There is no reason why the other predicted curves would 
not be seen in the experiment and hence the breakdown 
scenario is inconsistent with the data. In particular, al- 
though it could be argued that the exact value of B/Bq 
is uncertain and hence there is uncertainty in the rela- 
tive size of the BD2 and BDl/3 orbits, there is no good 
reason why the breakdown orbits at the top of the zone 
(BD2inax) should be seen and not that at the center of 
the zone (BD2niin)- The breakdown gap is not strongly 
kz dependent. 

Finally, we address estimates of the size of the break- 
down field. In Ref. {l| it is suggested that Bq is very 
small. Although even that does not resolve the problem 
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FIG. 3: Experimental dHvA data along with the ARPES 
bands (calculated in Ref. 1|), showing the non-breakdown 
orbits (dash lines) labelled as in Ref. [3] and the various pre- 
dicted breakdown orbits. The greek letters, e, 7 and /3 label 
the three experimentally observed frequency branches. 



of reconciling the dHvA and ARPES data, we would like 
to comment further on this. Bq can be estimated from 
the band-structure data using the following formula due 



to Chambers [11 j 
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where kg is the gap in fc-space, and a and b arc found from 
fits to the band dispersions close to the gap, ei — ep + 
fivp^{kx + ^kg+^aky) and 62 = eF+fiVF2{kx--^kg — ^bky). 
By fitting our band-structure calculations (see Fig. to 
this expression we get (at kz = 0.5c*) kg = 0.030A~^, 
a = 9.3 A and b = 8.6 A. This gives Bq = 125 T and 
hence a damping factor of approximately exp(— 500/i?) 
for the elliptic orbit (at k^ = we find Bq — 156 T). 
Clearly, this is much larger than our experimental fields 
explaining why we do not observe the magnetic break- 
down orbits. The measured damping factor for the e or- 
bit in our sample 0] is approximately exp(— 160/i3). This 
is consistent with impurity scattering alone (mean-free- 
path=400 A) , leaving no room for any significant contri- 
bution from magnetic breakdown. 

It is well known that correlation effects (which are rela- 
tively strong in LiFeAs where m* /rrib — 5 for some orbits) 
enhance SO and hence our estimated breakdown field is 
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FIG. 4: Fermi surface of the electron sheets in LiFeAs close 
to the point where magnetic breakdown is likely to occur (at 
kz = 0.5c*). The solid lines are fits to the quadratic expres- 
sion given in the text. 



likely too low. In the 122 compound it was found that the 
DFT calculations underestimate the magnetocrystalline 
anisotropy (spin gap) by a factor of two [l^ to three 
121 (meaning the SO coupling is underestimated). Also 
our analysis of LiFeP, where we have a nearly complete 
experimental picture, indicates that our calculations un- 
derestimate SO in the hole bands, and by implication, in 
the electron bands. In fact it is consistent with the fact 
that the ultra-small hole orbit 1 is not observed in LiFeAs 
(as opposed to LiFeP), if SO is stronger in LiFeAs. 

The ARPES spectra presented in the Comment fll| do 
not have sufficient resolution to determine this gap. Dif- 
ferent ARPES curves from the two electron bands are just 
about discernable at their maximum separation (124°) 
but as the SO gap in our calculations is 7.5 times 
smaller than this it would be impossible to see with the 
resolution shown in the figure. In addition, 'matrix' ele- 
ments effects mean that the second electron band is not 
observed for any angle less than 90°. It is unclear how the 
'matrix' elements evolve with angle so the SO gap may 
be unobservable even if there were sufficient resolution. 

To conclude, we would be very happy if our bulk dHvA 
were consistent with ARPES measurements, and in fact 
we had been looking for such an interpretation. Unfor- 
tunately, this does not appear to be borne out by the 
data. 
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